Spin dependent observable effect for free particles using the arrival time distribution 
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The mean arrival time of free particles is computed using the quantum probability current. This 
is uniquely determined in the non-relativistic limit of Dirac equation, although the Schroedinger 
probability current has an inherent non-uniqueness. Since the Dirac probability current involves 
a spin-dependent term, an arrival time distribution based on the probability current shows an 
observable spin-dependent effect, even for free particles. This arises essentially from relativistic 
quantum dynamics, but persists even in the non-relativistic regime. 
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I. INTRODUCTION 



The treatment of time in quantum mechanics is a much 
debated question. A testimony to this is the proliferation 
of recent papers [1] on the problems of tuneling time, de- 
cay time, dwell time and the arrival time. In this paper 
we are specifically concerned with the issue of arrival time 
[2]- 

In classical mechanics, a particle follows a definite tra- 
jectory; hence the time at which a particle reaches a given 
location is a well defined concept. On the other hand, 
in standard quantum mechanics, the meaning of arrival 
time is rather problematic. Indeed, there exists an ex- 
tensive literature on the treatment of arrival time dis- 
tribution in quantum mechanics [3]. A straightforward 
procedure would be to try to construct a self-adjoint op- 
erator for the arrival time in quantum mechanics which is 
conjugate to the Hamiltonian, but then it is seen that the 
operator does not have a basis of orthogonal eigenstates 
[4]- 

Using the Born interpretation, ti)\ 2 , |^(x, i 2 ) | 2 - - - 
give the position probability distributions at different 
instants t\, t^.... Now, the question posed is that if 
we fix the positions at x=Xi, X2..., can the functions 
|t/>(Xi, t)\ 2 , | , 0(X2, t)| 2 ...give the time probability distri- 
butions at different positions Xi, X 2 ...? Note that if at 
any instant t = tj, /_ °° \ijj(x.,t = ti)\ 2 d 3 x = 1, the prob- 
ability of finding the particle anywhere at that instant is 
unity. But if we fix the position at, say, x = Xi and t is 
varied, the value of the integral J* °° \ip(x = Xi,i)| 2 dt ^ 
1. In this case what may be pictured is that at a given 
point, say, Xi the relevant probability changes with time 
and this change of probability is governed by the follow- 
ing continuity equation which suggests a "flow of proba- 
bility" 



-|V(x,i)| 2 + V.J(x,i) =0 (1) 

where J(x, i)=^(V>VV>* — ip*Vip) is the probability cur- 
rent density. 

Different approaches for analysing the problem of ar- 
rival time distribution have been suggested using the 
path integrals and positive-operator- valued measures [5] . 
Delgado and Muga [6] proposed an interesting approach 
by constructing a self-adjoint operator having dimensions 
of time which is relevant to the arrival time distribution, 
but then its conjugate Hamiltonian has an unbounded 
spectrum. The implications of this approach have been 
studied in detail by Delgado [7]. 

In this paper we take recourse to the definition of ar- 
rival time distribution in terms of the quantum probabil- 
ity current density J(x = X, t). Interpreting the equation 
of continuity in terms of the physical probability flow, the 
Born interpretation for the modulus of the wave function 
and its time derivative seems to imply that the mean ar- 
rival time of the particles reaching a detector located at 
X is given by 

_ = f~\J(x = X,t)\tdt 

T / oo |J(x = X,t)| < ft U 

Here, it should be emphasized that the definition of the 
mean arrival time used in Eq.(2) is not a unique result 
within standard quantum mechanics. We also note that 
J(x, t) can be negative, hence one needs to take the mod- 
ulus sign in order to use the above definition. However, 
the Bohmian model of quantum mechanics in terms of 
the causal trajectories of individual particles leads to the 
above expression for the mean arrival time [8] in a rather 
conceptually elegant way. 

The quantum probability current interpreted as the 
streamlines of a conserved flux has been used for study- 
ing the tunneling times of Dirac electrons [9] . However, it 
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is easily seen that in non-relativistic quantum mechanics 
the form of the probability current density is not unique, 
a point which has already been discussed by a number 
of authors [10-12]. If we replace J by J' in Eq.(l) where 
J'=J+(5J, with V.<5J=0, J' satisfies the same probability 
conservation given by Eq.(l). Then this new current den- 
sity J' can lead to a different distribution function for the 
arrival time [12]. Hence the question arises how one can 
uniquely fix the arrival time distribution via the quan- 
tum probability current in the regime of non-relativistic 
quantum mechanics? 

In order to address the above question, we take a vital 
clue from the interesting result Holland [13] had shown in 
the context of analysing the uniqueness of the Bohmian 
model of quantum mechanics, viz. that the Dirac equa- 
tion implies a unique expression for the probability cur- 
rent density for spin 1/2 particles in the non-relativistic 
regime. In Section II we highlight the feature that the 
uniqueness of the probability current density is a generic 
consequence of any relativistic equation of quantum dy- 
namics. In Section III, the particular case of spin depen- 
dent probability current density as derived from the Dirac 
equation is discussed in detail. Subsequently, in Section 
IV, using the non-relativistic limit of the Dirac current 
density, we compute the effect of spin on the arrival time 
distribution of free particles for an initial Gaussian wave 
packet. Such a line of investigation has not been explored 
sufficiently; to the best of our knowledge, it seems that 
only Leavens [14] has studied this specifically in terms of 
the Bohmian causal model of spin-1/2 particles. 



Now, an interesting point is that this uniqueness will 
be preserved in the non-relativistic regime. Hence, given 
any relativistic wave equation, one can calculate the 
unique form of the current which can be applied in the 
non-relativistic regime. Thus using the (normalized) 
modulus of the probability current density as the arrival 
time distribution, if one calculates the mean arrival time, 
it can be used to empirically test any consistent relativis- 
tic wave equation such as the relativistic Kemmer equa- 
tion [16] for the massive spin and spin 1 bosons. Of 
late, the unique form of the probability current density 
expressions has been derived in the non-relativistic limit 
of the relativistic Kemmer equation for spin-0 and spin-1 
particles [17]. This general scheme for testing relativis- 
tic quantum wave equation in terms of the arrival time 
distribution is not contingent on any specific form of the 
relativistic wave equation. However, in the following de- 
tailed study we specifically use the Dirac equation for 
spin-1/2 particles. 



III. SPIN DEPENDENT EFFECT ON THE 
ARRIVAL TIME DISTRIBUTION USING DIRAC 
EQUATION 

The Dirac equation for a free particle is 



.^dib (he d . n\ , 
ih— = — cti- h f3m c \ip 



(3) 



II. UNIQUENESS OF THE PROBABILITY 
CURRENT DENSITY FOR ANY RELATIVISTIC 
WAVE EQUATION 

The current obtained from any consistent relativistic 
quantum wave equation will have to satisfy the covariant 
form of the continuity equation d^j^ = 0, where the ze- 
roth component of will be associated with the proba- 
bility density. Now, let us replace j M by f" which should 
again be conserved, i.e. d^j 1 " = 0, where = j^+a^, a M 
being an arbitrary 4-vector. But then the zeroth com- 
ponent j will have to reproduce the same probabil- 
ity density j , and hence o°=0. This current as seen 
from another Lorentz frame is j 11 ' =j IJ "-\-a^' . Then in this 
frame j° =j°+a° , and again from the previous argu- 
ment a =0. But we know that the only 4-vector whose 
fourth component vanishes in all frames is the null vector. 
Hence a M =0. Thus, for any consistent relativistic quan- 
tum wave equation, satisfying the covariant form of the 
continuity equation, the relativistic current is uniquely 
fixed. Unique expressions for the conserved currents have 
been explicitly derived by Holland [15] for the Dirac equa- 
tion, the Klein-Gordon equation, and also for the coupled 
Maxwell-Dirac equations. 



Cti = 



aA I 

a, P I -I 



ip is a four component column matrix and Oi are the 
Pauli matrices. Choosing a representation where 0i and 
■02 are two component spinors, one gets two coupled 
equations 



di>i dip2 im c 2 , 

~dt = ~ Ca ^ ~ 



d^2 dipi imoc 2 

-dt=- Ca ^ + —^ 



(4) 



(5) 



Combining Eqs.(4) and (5) one gets 

_( V , 1 ^ 1 ) = _ c ^ 1 V__ c __ (JiV , 1 ( 6 ) 

For positive energies, one can take 02 oc exp(—iEt/h). In 
the non-relativistic limit, E is the rest mass energy and 
E + m c 2 = (m + m )c 2 = 2m c 2 . Thus using Eq.(5) one 
can write 



1p2 



ihc dipi ih dtp 



a i^Zi=-^ZrZ a i^Zi ( 7 ) 



(E + m c 2 ) dx l 2m c dx 
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Putting this value of ip2 in Eq.(6), one gets 

V.J = 



dp 
dt 



(8) 



where J is the Dirac current in the non-relativistic limit 
that can be decomposed into two terms as was shown by 
Holland [13,15], as 



J = 



ih 
2~^i 



ih 



V'l^VV'i) - (V^i f )Vi 



^iV(ct.V)^i - (Vipi 1 .a)aipi (9) 

n 



+■ —V x (V»iWi) 
2m 



and p = ipi^ipi. ipi is a two component spinor which 
can be written for a particle in a spin eigenstate as 



1p! = 1p(x,t)x = 



R(x,t) exp 



iS(x, t) 
h 



X 



(10) 



Here tjj(x, t) is the Schroedinger wavefunction and x 1S 
a spin eigenstate. Putting this form of ipi in the expres- 
sion for current in Eq.(9) one gets 



J = — pVS+ — CVp x s) 
m m 



with 



s = (h/2) X f a X 



(11) 



R 2 



The first term (Jj) in Eq.(ll) is independent of spin, 
while the second term (J s ) is the contribution of the spin 
of a free particle to the unique conserved vector current in 
the non-relativistic limit. It is then clear that the mean 
arrival time given by Eq.(2) can be computed by using 
the unique expression for J in Eq.(ll). Thus we can ob- 
tain a spin-dependent contribution in the expression for 
the mean time of arrival for free particles, which could 
be an experimentally measurable quantity. On the other 
hand, by ignoring the spin-dependent term one would 
obtain the mean arrival time given by 



n = 



Jq 00 \3i\tdt 



(12) 



In the following section IV we study the situations where 
the difference between the actual magnitudes of f and fj 
is significant, thereby enhancing the feasibility of exper- 
imentally testing the specific spin-dependent effect. 



IV. THE COMPUTED EFFECTS ON THE 
ARRIVAL TIME DISTRIBUTION 

We consider a freely evolving Gaussian wave packet 
in the two separate cases (A and B) corresponding to 



an initially symmetric and an asymmetric wave packet 
respectively. 

Case A: Symmetric wave packet 

Let us consider a Gaussian wave packet for a free spin 
1/2 particle of mass m centered at the point x = 0, 
y = 0,and z = 0. We choose the spin to be directed 
along the z-axis, i.e., (s = jz). 



V(x,i = 0) 



1 



(2tr7 2 ) 3 /' 



-exp(ik.x)exp 



The time evolved wave function can be written as 

~iS(-x,t)~ 



ip(x, t) = R(x, i)exp 



where 



i?(x, t) = (2Tra 2 ) 3/4 exp 



(x - ut) 2 



4cr 2 



and 



ht 



S(x,t) = -^tan" 1 ( n 
K ' 2 \2mao 2 



(13) 



(14) 



(15) 



(16) 



1 , (x-ut) 2 h 2 t 

+mu.(x - -ut) H 5-5 — 

2 StticqC 

with (u = hk/m) the initial group velocity taken along 
the x-axis, and 



a = (T 



2+2 



1 1/2 



4m 2 <4 



(17) 



The total current density can be calculated using Eq.(ll) 
to be (we set m = 1, h=l) 



(x — ut)t 
u + — „ / |x 



yt 



4a V2 



y + 



zt 



+ p 



2a 2 



x 



(x — ut) „ 
2a 2 Y 



(18) 



where the contribution of spin is contained in the second 
term only. 

We can now compute f and fi numerically by sub- 
stituting Eq.(18) in Eqs.(2) and (12) respectively. It is 
instructive to examine the behaviour of the contribution 
of spin-dependent term towards the mean arrival time. 
For this purpose, we define a quantity 



_ I °°\J s \tdt 

Ts Cwt 



(19) 



We first compute f s for a range of the initial velocity u in 
units of m=l, and U=l. We find that the spin of a free 
particle contributes towards altering its mean arival time 
for a wide range of initial velocities. This feature holds 
generally, except for very small magnitudes of velocity 
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where the spin-dependent contribution may be negligi- 
ble depending on the location of the detector vis-a-vis 
the direction of the initial group velocity u. This fea- 
ture is shown in Figure. 1 where we plot the variation of 
f s with u. The initial wave packet is peaked at the ori- 
gin with (To = 0.01. The detector position is chosen at 
(x = l,y = 1, z = 1). We also find that the difference of 
magnitude between f and Tj can be increased by choos- 
ing asymmetric detector positions as well as asymmetric 
spread for the initial wave packet, an example of which 
we highlight below. 



0.00024 



0.00022 - 



0.0002 



0.00018 



0.00016 



0.00014 



where a = (a 2 + it) 1 / 2 : (3 = (b 2 + it) 1 ' 2 - 7 = (c 2 + it) 1 ' 2 . 
Writing the wave function as 



ip(x, y, z, t) = R(x, y, z, t)exp 
one obtains 



iS(x,y,z,t) 



R(x,y,z,t) 



aW^ 1/4 



l 



(22) 



0.00012 




exp 
exp 



a 2 (x + x\ - ktf- 
2(a 4 + t 2 ) 



exp 



(p2 + g 2)l/4 

b 2 y 2 



2(fo 4 + i 2 ) 



2(c 4 + t 2 ) 



(23) 



and 



S(x,y,z,t) = hkx - - ^tan \q/p) 



ht(x + X! - ktf 
2(a 4 + t 2 ) 



+ 



My 2 



2(b 4 + t 2 ) 2(c 4 + i 2 ) 



(24) 



with 



p = a 2 b 2 c 2 - a 2 t 2 - b 2 t 2 - c 2 t 2 



250 



q = a 2 b 2 t + a £ c l t + b 2 c l t - t 



2„2, , >2 2. 



FIG. 1. The spin-dependent contribution to the mean ar- 
rival time computed at the point x=l, y=l, z=l is plot- 
ted against the initial group velocity of the packet along the 
x-axis. 



Case B: Asymmetric wave packet 

We consider an initial free particle wave packet in three 
dimensions which is centered at the point x = —Xi, 
y = 0,z = 0. 



(25) 



Considering again a spin-1/2 particle with spin di- 
rected along z-axis (s = ^z), the total current density 
defined in Eq.(ll) is given by (in units of h = 1 = m) 



J = p 

+P 



(x + x\ — ut)t s 
U+ (a 4 + t 2 ) iX 



yt 



zt 



(b 4 + t 2 y (c 4 + t 2 ) 



b 2 y „ a 2 (x + x\ - ut) „ 
:x+ , A . — y 



(b 4 + t 2 Y 



(a 4 + t 2 ) 



(26) 



where the second term represents the spin-dependent 
contribution to the current. 



ip(x,y,z,t= 0) 



1 



ir 3 a 2 b 2 c 2 



1/4 



exp(ikx) 





{X + Xx) 2 ' 




[ (y) 2 ' 




\ 


exp 


2a 2 


exp 


2b 2 


exp 


2c 2 



(20) 



where a, b, c are positive constants. (Such a form for the 
wave function was considered by Finkelstein [12] in the 
context of arrival time distributions.) The particle is 
given an initial velocity in the x direction represented 
by u — —. The time evolved wave function is given by 



ip(x,y,z,t) = 
(x + X\ 



a 2 b 2 c 2 ^ 1/4 



exp 



ktf 



exp[i(kx - k 2 t/2)} 



0.865 



0.855 - 



0.85 



0.845 



0.84 



0.835 



2a 2 



exp 



2/3 2 



exp 



2 7 2 



(21) 
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FIG. 2. The mean arrival times f (upper curve) and f< 
(lower curve) computed at the point x = l,y = 2,z = 1 are 
plotted against the initial group velocity of the packet along 
the x-axis. and a = 0.001, b = 0.4, c = 0.01, xi = 0. 



We compute numerically the arrival times f and f». 
Figure. 2 shows the variation of f and fi with the initial 
group velocities (u) of the wave packet. Here we choose 
the parameters as x\ = 0,a = 0.001,6 = 0.4, c = 0.01. 
Accordingly the mean arrival time is computed at the 
position x = 1.0, y = 2.0, z = 1.0. One sees that the dif- 
ference in the magnitudes of f and fj can suitably be en- 
hanced by a judicious choice of asymmetric initial spreads 
and detector positions. 

V. CONCLUDING REMARKS 

Let us now summarise the salient features of our 
scheme. For measuring the spin of a particle, it is usu- 
ally subjected to an external field, like in a Stern-Gerlach 
apparatus. But the scheme we have suggested would 
enable to detect, a spin-dependent effect without using 
any external field. Such an observable effect thus high- 
lights the feature that the spin of a particle is an intrinsic 
property and is not contingent on the presence of an ex- 
ternal field. As demonstrated in this paper, the spin- 
dependent term in the Dirac probability current density 
contributes significantly to the computed mean arrival 
time for a range of suitably chosen parameters of the 
Gaussian wave packet. Thus if the arrival time distri- 
bution can be measured, this predicted spin-dependent 
effect would be empirically verifiable. 

Another way of perceiving the significance of such an 
effect is as follows. The dynamical properties of free par- 
ticles like position, momentum, and energy can of course 
be measured. However, one cannot usually measure the 
static or innate particle properties such as charge with- 
out using any external field. Nevertheless, the scheme 
we have proposed shows that the magnitude of total spin 
can be measured without subjecting the particle to an 
external field. 

Another interesting implication of the measurement of 
the spin dependent arrival times for free particles could 
be to view this as implying a fundamental difference be- 
tween the magnitude of total spin of a particle and its 
other static properties such as mass and charge in the fol- 
lowing sense. The property of spin seems to be crucially 
contingent on the fundamentally rclativistic nature of the 
dynamical equation of the wave function of the particle 
so that the wave function is essentially 4-componcnt, or 
2-component in the non-relativistic limit. 

Here we would like to stress that the spin-dependent 
term which contributes significantly to the arrival time 



distribution we have computed in the nonrelativistic 
regime originates from the relativistic Dirac equation and 
hence this provides a rather rare example of an empiri- 
cally detectable manifestation of a relativistic dynamical 
equation in the nonrelativistic regime, an effect which 
cannot be derived uniquely from the Schroedinger dy- 
namics. 

A further line of investigation as an offshoot of this pa- 
per could be to explore the possibilities of using the rela- 
tivistic quantum mechanical wave equations of particles 
with spins other than spin 1/2 (such as using the Kemmer 
equation [16,17] for spin and spin 1 particles) in order to 
compute the spin-dependent terms in the probability cur- 
rent densities and their effects on the arrival time distri- 
bution. Such a study seems worthwhile because then the 
arrival time distribution may provide a means of checking 
the validity of the various suggested relativistic quantum 
mechanical equations which have otherwise eluded any 
empirical verification. 
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